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Abstract 

An Abelian sandpile model is considered on the Husimi lattice of square pla- 
quettes. Exact expressions for the distribution of height probabilities in the Self- 
Organized Critical state are derived. The two-point correlation function for the 
sites deep inside the Husimi lattice is calculated exactly. 
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1 Introduction. 



In recent years, there has been considerable interest in different dynamical models which 
can evolve without any tuning of parameters to the Self- Organized Critical (SOC) state. 

The concept of Self- Organized Criticality has been introduced by Bak, Tang, and 
Wiesenfeld |] , through simple cellular automaton models known as sandpiles, to explain 
the temporal and spatial scaling in dynamical dissipative systems. 

Later on, Dhar and Majumdar, in a number of articles [2-4], have studied the so- 
called Abelian sandpile models (ASM's) for they show a nontrivial analytically tractable 
example of the SOC behavior. For these models some ensemble-average quantities have 
been calculated on the square lattice [2,4-7]. 

In this article, we consider ASM on the Husimi lattice of square plaquettes. One of 
the remarkable features of the Husimi- or Bethe-like lattices is the exact solvability of 
different spin, gauge and dynamical models defined on them [8-14]. One might also hope, 
by examining the form of exact solutions for these lattices which allow only a limited 
number of closed configurations, to predict a general behavior for the lattice models of 
principal interest. 

It is well known that the exact solutions obtained for different spin models on the 
Bethe lattice are equivalent to the Bethe-Peierls approximation and they give a more pre- 
cise description of the critical behavior of these models than the mean-field approximation 
technique ||15||. In turn, the Husimi lattice can be considered as a next step in this list of 



approximations and, as it has been shown by Monroe |L0| for the two- and three-site in- 
teracting spin systems, this approximation improves results obtained on the Bethe lattice. 
In section 6 we bring a table where we compare the distribution of height probabilities in 
the SOC state of ASM calculated on different lattices and again we see that the results 
obtained on the Husimi lattice are in good agreement with the known exact results on 
the square lattice @, || . This can be explained by the presence of elementary loops in the 
Husimi lattice that reproduce the local structure of usual lattices more precisely than the 
Bethe lattice. Hence, for most of flat and three-dimensional lattices we can find suitable 
Husimi-like lattices and examine different models on them. 

The outline of the article is as follows: in the next section, we define the lattice 
and ASM on it. In section 3, we find recursion relations for the numbers of allowed 
configurations in the SOC state. In section 4, we exactly compute the distribution of 
height probabilities. In section 5, we calculate the two-point correlation function for the 
sites deep inside the lattice. In section 6, we bring some concluding remarks. 



2 Lattice and Model. 

A pure Husimi tree of square plaquettes |1| can be constructed recurrently. As a basic 
building block, we will take an elementary plaquette (Fig. la). This basic block will 
be called the first- generation branch. To construct the second-generation branch, we 
will attach a single basic building block at each free site of the first-generation branch 
except the base site (root) (Fig. lb). Continuing this process we develop higher-generation 
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branches. Then, at the final step, we will take four nth-generation branches and connect 
their base sites by the elementary plaquette. As a result, we will get the graph with the 
coordination number q = 4. 

Let as define ASM on this connected graph of iV sites as follows: to each site % 
(1 < % < N) we associate an integer Zi (1 < z^ < 4) which is the height of a column of 
sand grains. The evolution of the system is specified by two rules: 

i. Addition of a sand grain at a randomly chosen site i increases Zi by 1. 

ii. The site i topples if the height Z{ exceeds the critical value z c = 4 and sand grains 
drop on the nearest neighbors. 

The number of surface sites of the Husimi tree is comparable with the interior ones. 
Hence, the calculation of the thermodynamic limit of the bulk properties requires special 
care. In our work, we define the height distribution of sand grains and the two-point 
correlation function for the sites deep inside the tree. Using these interior sites one can 
construct an infinite lattice, as they have the same features. Therefore, we will consider 
the problem on the Husimi lattice rather than on the Husimi tree. 

Any configuration {z^} on the Husimi tree in which 1 < < 4 is a stable configuration 
under the toppling rule. These configurations can be divided into two class: allowed and 
forbidden configurations M . 

In the SOC state, only allowed configurations have a nonzero probabiliy. Any subcon- 
figuration of heights F on a finite connected set of sites is forbidden if 

Zi < qu V i G F , 

where qi is a coordination number of a site i in the given subconfiguration F B . 

In turn, we can divide the allowed subconfigurations on an nth-generation branch of 
the Husimi tree into three nonoverlapping classes: weakly allowed of the type 1 (Wi), 
weakly allowed of the type 2 (W 2 ) and strongly allowed (S) subconfigurations. 

Consider an allowed subconfiguration C on the nth-generation branch G n with a root 
a (Fig. 2a). The coordination number of the root a is q = 2. Adding a vertex b to the G n , 
one defines a subgraph G' = G n U b. If the subconfiguration C = C U b with z b = 1 on G' 
is forbidden, then C is called the weakly allowed subconfiguration of type 1 (Wi). Thus, 
W\ can be locked by one bond, after which it becomes forbidden. 

Now add two vertices b and d to G n and consider a subconfiguration C" = C U b U d 
on G" = G n U b U d (Fig.2b). If the subconfiguration C" = C U b U d with z b = 1 and 
Zd = 1 on G" is forbidden, then C is called the weakly allowed subconfiguration of type 2 
(W 2 ). 

Any allowed subconfigurations defined on the nth-generation branches that cannot be 
locked by one bond or by two bonds form a strongly allowed (S) class. 

It is important to note that any subconfiguration of the W\ type is also of the W 2 
type. To obtain the nonoverlapping classes, we always check the subconfiguration first to 
belong to the W\ type and only then to the W 2 type. 
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3 Recursion Relations. 

In this section, we endeavor to find recursion relations for the numbers of allowed con- 
figurations on the nth-generation branch G n of the pure Husimi tree. Consider now G n 
with a root vertex a (Fig.3) that consists of three (n — l)th-generation branches G^li, 
G^-i and G^-\ with roots a l5 a 2 and a 3 , respectively. Let N Wl (G n ,i), N W2 (G n ,i) and 
Ns(G n ,i) be the numbers of distinct W±, W 2 and S type subconfigurations on G n with a 
given height z a = % at the root a. 
Let us also introduce 

4 

N Wl (G n )=Y,N Wl (Gn,i), (3.1) 
i=i 

N W2 (G n ) = j2 N w*(Gn,i), (3.2) 
i=i 

N s (G n ) = j2N s (G n ,z). (3.3) 
i=i 

These numbers can be expressed in terms of the numbers of allowed subconfigurations 
on the three (n — l)th-generation branches G^li, G^li and G^li- 

N Wl (G n ) = N S N S N S + N Wl N s N s + N W2 N S N S + N s N Wl N s + N W2 N Wl N s + 

N S N W2 N s + N Wl N W2 N s + N W2 N W2 N s + N s N s N Wl + N Wl N s N Wl + (3.4) 
N W2 N s N Wl + N s N W2 N Wl + N W2 N W2 N Wl + N S N S N W2 + N Wl N s N w , 2 + 
N W2 N S N W2 + N s N Wl N W2 + N W2 N Wl N W2 + N S N W2 N W2 + N Wl N w , 2 N w , 2 + N W2 N W2 N m , 

N W2 (G n ) = N Wl (G n ) , (3.5) 



N s (G n ) = 2N S N S N S + N Wl N s N s + 2N W2 N S N S + 2N s N Wl N s + N W2 N Wl N s + 

2N S N W2 N S + N Wl N W2 N s + 2N W2 N W2 N S + N s N s N Wl + N W2 N s N Wl + (3.6) 
N s N W2 N Wl + 2N S N S N W2 + N Wl N s N W2 + 2N W2 N S N W2 + N s N Wl N W2 + 2N S N W2 N W2 , 

where the first factor in each term of the sum corresponds to the G n l \ branch, the second 

(2) (3) 

one to and the third one to G n Li- 

The fact that the numbers of the W\ type subconfigurations and the W 2 type ones are 
equal to each other is seen from the straightforward calculation. 
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Let us introduce 

* = < 3 - 7 > 

where N w = N Wl = N W2 . 

If we consider graphs G^li and G„2i to be isomorphic, then iV(G„_i) = 

^(Gn-i) = ^(Ci-i) and from Q -(13.71) one obtains the following recursion relation 



l+4X(G ra _ 1 ) + 2X 2 (q t _ 1 ) 
X(Gn) " 2[1 + 3X(G B _ 1 )] • (3 - 8) 
The iterative sequence {X{G n )} starts from the seed X(G ) = ~ and converges to the 
stable point X* = 1+ 4 V ^ that characterizes, in the thermodynamic limit, the ratio of the 
W\ type or the W2 type configurations to the strongly allowed ones. 



4 Distribution of Height Probabilities. 

One of the main characteristics that describes the SOC state is the probability P(i) of 
having the height z = i at a given site: 

= (4.D 

jVtotal 

where N(i) is the number of allowed configurations with a given value z — i (1 < i < 4) 
and iVtotai = Si=i -^(«) is the total number of allowed configurations on the Husimi lattice. 

Consider now a randomly chosen site O deep inside the Husimi tree (Fig. 4). The 
number N(i) can be expressed via the numbers of allowed configurations on the six nth- 
generation branches Gffl, a = 1, ... ,6. 

If i = 1, then each allowed configuration on the branches G^\ G^\ G$ and G^ 
cannot be of the W\ type. It is also evident that three W2 type configurations cannot 
occur on the neighboring branches G^\ G^\ G^\ or G^\ G^\ Gffl, and so on. 

Excluding all such forbidden subconfigurations one finds the following expression for 
the number Nil) of allowed configurations with i = 1 at O for isomorphic branches G^\ 
a = 1, . . . , 6: 

6 

N(l) = [1 + 8X + 22X 2 + 24X 3 + 9X 4 ] J| X s (G (a) ) . (4.2) 

a=l 

Arguing similarly, one can get 



N{2) = [1 + 12X + 50X 2 + 84X 3 + 45X 4 ] f[ N s (G {a) ) , (4.3) 

a=l 

6 

X(3) = [1 + 12X + 56X 2 + 124X 3 + 119X 4 + 24X 5 ] X s (G (a) ), (4.4) 

a=l 
6 

N(4) = [1 + 12X + 56X 2 + 128X 3 + 147X 4 + 72X 5 ] J[ N s (G {a) ) . (4.5) 



0=1 
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For the sites far from the surface in the thermodynamic limit (n 
X = i±^I. Thus, from( fO)-§j|) we § et 



oo) we have 



P(D 



5(>/5 — 2) 



16 



P(2) 



10-3^5 
16 



P(3) 



3(4- >/§) 



16 



4 + 
16 



(4.6) 



5 Two-Point Correlation Function. 



The two-point correlation function may be defined as the probability P n (i,j) of a stable 
configuration in which two sites separated by the distance n have heights i and j. To 
obtain P n (i,j) explicitly, we consider two sites A\ and A n+ i deep inside the Husimi lattice 
(Fig.5). The left-hand side of the lattice beginning from the vertex Ak, k — 1, . . . , n, will 
be denoted as branch or subtree Gk with the root A k . This branch, in turn, consists of 
the subbranches Gk-i, U% , U^li- 

Following Dhar and Majumdar ||, we can solve the problem by the transfer matrix 
technique based on the fractal structure of the lattice. 

The number of allowed configurations N(i,j) on the Husimi lattice with fixed heights 
i at Ax and j at A n+1 is expressed via N a (G n ), N a (U®), N a (U^ ] +l ), N a (U n+2 ), N a (uS+ 2 ) 
and N a (Un+i), where a = W\, W 2 , S. The bars above N a (G n ) denote the constraint that 
the height at A\ is fixed at i. 

The numbers N s (G n ), N Wl (G n ) and N W2 (G n ), in turn, can be expressed in terms of 
the numbers of allowed configurations on the branches U^lx and G ra _i in the matrix 
form 

n-l 



( N Wl (G r 
Nw 2 (G r 
V N s (G n 



N S {US») • N s (U t 



(6) ■ 
n— 1, 



A 



Nw 2 (Gn-l) 
V N S (Gn-l) 



where 



A 



/ 1 + 3X + X 2 
1 + 3X + X 2 
V 1 + 2X 



1 + AX + 3X 2 
1 + AX + 3X 2 
2 + QX 



1 + AX + 3X 2 

1 + AX + 3X 2 

2 + 7X + AX 2 



In the thermodynamic limit for the sites deep inside the lattice we have X 
Hence, one can get 



(5.1; 



(5.2) 



/ N Wl (Gn 

Nw 2 {G n 
V Ns(G n ) 



HNs{ujl») • N 8 {u!t>) - A 



n-l 



/ N Wl {G u i) \ 
N W2 (Gx,i) 



(5.3) 



where 



/ N Wl (Gx,i) 
N W2 {Gi,i) 
V NsiGxJ) 



( NwAGi) 
Nw 2 {Gx) 
\ N s (Gi) 



and the product runs over all subbranches of the subtree G n . 



(5.4) 
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To obtain the two-point correlation function, we ought to divide N(i,j) by the total 
number of allowed configurations of the full lattice in the SOC state 



PndJ) = §^, (5.5) 



where iV total is given by 



iVtoui = 521 + 2 233 ^ A r 5(G g iVg (^ (5.6) 

with ^ 

N s (G n )=J2Ns(G n ). (5.7) 
i=i 

Then, substituting N s (G n ) from ( [5 . 3| ) one can find for n > 1 

iVtotal = 521 ^ 1165 (18 VE + 40 + (3 + v^XAa - A 2 )) X^ 1 J[ N S (U) , (5.8) 
where 

Ax = 0, (5.9) 



21 + 9v / 5>V /470 + 210v / 5 
A 2 ,3 = \ , (5.10) 

are eigenvalues of the matrix A and the product runs over all subbranches which surround 
the path from A\ to Ai+i- 

Thus, after a rather tedious calculation we obtain an exact expression for P n (i,j) in 
the SOC state 

/A \ _(n_1) 

P n (iJ) = P(i)P(j)+ Pij Utj , n>l, (5.11) 

where numerical constants 

= 8140 + 3640 - (1282 V$ + 2870) (A 3 - A 2 ) , 

= 50660 + 22656 - (2762 V5 + 6174) (A 3 - A 2 ) , 

: 122760 + 54900 V$ - (3234 Vh + 7230) (A 3 - A 2 ) , 

: 171560 + 76724 - (3026 + 6766) (A 3 - A 2 ) , 

= 175612 + 78536 V$ - (— ==- + 14982)(A 3 - A 2 ) , 

V5 

= 333040 + 148940 - (9290 + 20774) (A 3 - A 2 ) , 
r- 49922 

: 425488 + 190284 \fh - ( — -=- + 22326) (A 3 - A 2 ) , 

V5 
76682 

= 521612 + 233272 - ( — 7=- + 34294) (A 3 - A 2 ) , 

V5 





Pll 


Pl2 


= P21 


Pl3 


= P31 


Pl4 


= P41 




P22 


P23 


= P32 


P24 


= P42 




P33 



p 3 4 = p 43 = 605724 + 270888 V$ 



,91674 



+ 40998) (A 3 - A 



2) , 



r ,115466 

p 44 = 662860 + 296440 Vb - ( + 51638)(A 3 - A 2 ) . 

V5 

For two nearest-neighbor sites A\ and A 2 separated by the distance n — 1 in the SOC 
state the constants p^- have other values because the matrix A is degenerate. Hence we 
get 



Px(3,3) 
P 1 (1,2)=P 1 (2,1) = 
P 1 (1,3) = P 1 (3,1) = 
P 1 (l,4) = Pi(4,l) = 

6 Conclusion. 



Pi(l,l) = 

119-47^5 
" 128 
5(13^5-29) 
32 

265 - 11A/5 
128 

235 - 103^ 
128 



Pi(2,2) 
A (4, 4) 



175 - 78^ 
16 

87-31^ 



128 



P x (2,3) = Pi(3,2) 
P 1 (2,4) = P 1 (4,2) 
P 1 (3,4) = P 1 (4,3) 



169^-369 

128 
171^5-371 

128 
81 - 29^ 
128 ' 



In this article, we have investigated the Abelian sandpile model on the Husimi lattice of 
square plaqettes. The distribution of height probabilities and the two-point correlation 
function in the SOC state have been calculated exactly. We have shown that correlations 
decay with distance as (A 3 /A 2 )~( n_1 ). 

At the end we want to compare the distribution of height probabilities obtained on 
different lattices with the same coordination number q = 4: 





square 0, || 


Bethe § 


Husimi of triangles 


14 




Husimi of squares 


P(l) 


0.07363 


0.07407 


0.07031 


0.07377 


P(2) 


0.1739 


0.22222 


0.16406 


0.20574 


P(3) 


0.3063 


0.33333 


0.33594 


0.33074 


P(4) 


0.4461 


0.37037 


0.42969 


0.38975 



From this table we see that the Husimi lattice is a good approximation for the regular 
lattices and the best one is achieved for P(l) as it can be considered as a local structural 
characteristic of the model ||. The next step of our investigations will be the calculation 
of dynamic characteristics of ASM, e.g., critical exponents of avalanches. The choice of 
the Husimi lattice gives us hope to find the relationship between the chaos and SOC state 
since some spin models formulated on this lattice show the chaotic behavior |T0], O . 
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base site 



base site 



Figure 1: (a) A first-generation branch consists of a single square plaquette. (b) A second- 
generation branch. 
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a 



Figure 3: The nth-generation branch G n with three nearest (n — l)th-generation branches 




Figure 4: A site O with height i is located deep inside the lattice and surrounded by the 
six nth-generation branches Gg\ Gg\ Gf \ G<$, and G%\ 




Un+2 




U n+2 



Figure 5: The path from the site A\ to the site A n+ i on the Husimi lattice goes through 
the points A 2 ,. ■ ■ , A n . The left-hand side of the lattice beginning from the vertex A k , 
k — 1, . . . , n, is denoted as a branch Gk with the rootA^. 
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